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(Dated: November 1, 2012) 

We suggest that the combined effect of screening, gluon-induced dissociation, collisional damping, 
and reduced feed-down explains most of the sequential suppression of T(nS) states that has been 
observed in PbPb relative to pp collisions at ^snn = 2.76 TeV. The suppression is thus a clear, albeit 
indirect, indication for the presence of a QGP. The T(1S) ground state suppression is essentially 
due to reduced feed-down, collisional damping and gluodissociation, whereas screening prevails for 
the suppression of the excited states. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Cj 



I. INTRODUCTION 

The suppression of quarkonium (QQ) states is one 
of the most promising probes for the properties of the 
quark-gluon plasma (QGP) that is generated in heavy- 
ion collisions at highly relativistic energies. In the QGP 
the confining potential is screened due to the interaction 
of the heavy QQ with medium partons and hence, char- 
monium and bottomium states successively melt [l[ at 
sufficiently high temperatures Tdi SS beyond the critical 
value T c m 170 MeV. 

However, additional processes such as gluon-induced 
dissociation, and collisional damping contribute to the 
suppression, and are effective in a temperature region 
where the T(nS) states - and in particular, the T(15) 
ground state - have not yet melted due to screening. 

Here we concentrate on such processes. It turns out 
that in particular for the T(15) ground state, bot- 
tomium dissociation is not just static screening, but 
mostly caused by other means - whereas the dissocia- 
tion of the excited states is essentially due to screening. 

Charmonium suppression has been studied since 1986 
in great detail both theoretically and experimen- 

tally at energies reached at the CERN SPS, BNL RHIC 
[j| , and the CERN LHC @, 0] ■ Bottomium suppression is 
expected to be a cleaner probe. The T(15) ground state 
with mass 9.46 GeV is strongly bound. It melts as the 
last QQ in the QGP (depending on the potential) only at 
about 4.10 T c Even at LHC energies the number of 
66-pairs in the QGP remains small such that statistical 
recombination is unimportant. 

T suppression in heavy-ion collisions has recently been 
observed for the first time both by the STAR experiment 
at RHIC [|, and by the CMS experiment at LHC flllflljj. 
Preliminary CMS data from the 2011 run (l2j have much 
better statistics such that the T(25) state can now be 
resolved individually in PbPb collisions at the LHC. 

In this work we suggest a three step model that consid- 
ers the T (IS, 25, 35) and XbO-P, 2P) states to obtain the 



suppression of the T(15, 25, 35) states at LHC energies, 
which is then compared to the experimental results. We 
successively calculate 

1. the bb wave functions, and decay widths for the 
three processes Debye screening, collisional damp- 
ing and gluodissociation [l| 



2. the suppression of the five states considered here 
within the expanding and cooling fireball 



3. the feed-down cascade, and the ensuing fraction of 
dimuon decays, T(n5) 
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Whereas gluodissociation below T c is not possible due 
to confinement, it does occur above T c where the color- 
octet state of a free quark and antiquark can propagate 
in the medium. Its significance increases substantially 
with the rising gluon density at LHC energies. 

In the midrapidity range \y\ < 2.4 where the CMS 
measurement [10h12| ] has been performed, the tempera- 
ture and hence, the thermal gluon density is high, and 
causes a rapid dissociation in particular of the T(25) 
and T(35) states, but also of the T(15) ground state. 
At larger rapidities up to the beam value of j/bcam = 7.99 
and correspondingly small scattering angles where the 
valence-quark density is high [l4| . nonthermal processes 
would be more important than in the midrapidity region 
that we are investigating here. 

The calculation of the bb wave functions for five T(n5) 
and \b(nP) states, and the associated widths Tdamp of 
these states due to collisional damping from a complex 
potential are considered in the following section. The cal- 
culation of the gluodissociation decay widths Tdiss for the 
same states is discussed in Sec. Ill, the time evolution of 
the fireball and subsequent decay cascade is considered 
in Sec. IV. The results are presented in Sec. V in com- 
parison with the available CMS data at LHC energies, 
and the conclusions are drawn in Sec. VI. 
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II. BOTTOMIUM WAVE FUNCTIONS AND 
COLLISIONAL DAMPING 

Due to the small relative velocity v <C c of the bottom 
quarks in the bound state, bb may be properly described 
by the potential NonRelativistic QCD (pNRQCD) ap- 
proach [HI, [l6| ■ The relevant terms in the pNRQCD ac- 
tion for the 66-pair read [see e.g. 
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with the singlet and octet fields S and O a , the ultra 
soft color electric field E a , the bottom mass m b — 4.89 
GeV, the strong coupling constant at the soft scale, a s s — 
a s {m b a s /2) = 0.48, and N c = 3, C F = {N 2 c -l)/ \2N C ) = 
4/3. 

This approach leads to a Schrodinger equation, with 
the coulombic, color-singlet potential V = —Cpa s s /r. 
For the treatment of bb in the QGP it is, however, appro- 
priate to make a calculation at finite temperature which 
yields for the short-range part of the potential, in the 
HTL approximation, a complex, screened, coulombic ex- 
pression Hi m. 

The long range part is parameterized as in [2 31 ] so that 
the full singlet potential reads 
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with a c tf = Aa s s /3, the Debye mass tod, number of flavors 
in the QGP Nf — 3, and the strong coupling constant 
evaluated at the HTL energy 2ttT, a T s = a s {2irT) < 0.50, 
respectively. The absolute values |(?ni(^)| of the resulting 
bb wave functions are shown in Fig. [TJ 

The Schrodinger equation is now solved for every bb 
state with the potential @ for T > T c up to the dissoci- 
ation temperature Tdi SS above which screening prevents 
bottomium formation and the Schrodinger equation has 
no bound states solutions. The dissociation temperatures 
with the above parameters are Tdiss — 668, 217 and 206 
MeV for the T{1S), T(2S) and x b (lP), respectively: The 
higher excited states are already dissolved for T >T C . 

The imaginary part of the potential causes a decay 
width Tdamp, monotonically increasing with temperature, 
which accounts for collisional damping by the plasma 
particles. Tdamp is displayed in Fig. [3] together with the 
decay width Tdiss for gluodissociation, which is consid- 
ered in the next section. 
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FIG. 1. (Color online) Radial wave functions of the T(15), 
Xb(lP) and T(2S) states (solid, dotted, dashed curves, re- 
spectively) calculated in the complex, screened potential ((2} 
for temperatures T = MeV (bottom) and 200 MeV (top) 
with effective coupling constant a e fl = (4/ 3) a* = 0.63, and 
string tension a — 0.192 GeV 2 . While the rms radius U (r 2 ) 
of the T(IS') ground state is almost insensitive to tempera- 
ture changes, it varies substantially with temperature for the 
Xt(lP) and T(25) states. 



III. GLUODISSOCIATION IN THE MEDIUM 

Due to the high gluon density reached at LHC energies 
in the mid-rapidity region, gluodissociation is a major 
process besides screening and collisional damping that 
leads to a suppression of T's at LHC. Hence we calculate 
the gluodissociation cross sections for the T(1S)-T(35), 
and XbO-P), Xb(2P) states for different lifetimes £qgp of 
the QGP. 

The leading-order dissociation cross section of the bb 
states through El absorption of a single gluon had been 
derived by Bhanot and Peskin (BP) [24(. From the pN- 
RQCD approach the gluodissociation cross section may 
be derived from the dipole interaction term in eq. ([1} de- 
scribing a singlet-octet transition of the 66-pair via emis- 
sion/absorption of an ultra soft gluon. From this starting 
point we can easily generalize the approach to include the 
effect of our modified potential ^ [25| , and obtain for a 
bottomium state (nl) 
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FIG. 2. (Color online) Gluodissociation cross sections 
fdiss ( n S) in mb (left scale) of the T(1S) and T(2S) states 
calculated using the screened complex potential for temper- 
atures T = 170 (solid curves) and 200 MeV (dotted curves) 
as functions of the gluon energy E g . The thermal gluon dis- 
tribution (right scale; solid for T = 170 MeV, dotted for 200 
MeV) is used to obtain the thermally averaged cross sections 
through integrations over the gluon momenta. 
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with the singlet and octet states \nlm), \ql'm!) and a" = 
a s (rriba 2 /2) ~ 0.59. The radial wave function h q i> of the 
states \ql'm'} is derived from the octet Hamiltonian with 
the potential V$ = +a e ff/(8r). 

We had originally derived the gluodissociation cross 
section in [l3| independently from the pNRQCD formu- 
lation in an approach that was based on a straightfor- 
ward extension of the Bhanot-Peskin formulation [24( to 
approximately account for the confining string contribu- 
tion m. 

For vanishing string tension and the corresponding val- 
ues of the binding energy E n i, a pure Coulomb IS wave 
function, and a simplification in the octet wave function, 
our expression reduces to the result in [24]]. Our full re- 
sult for the T(IS') gluodissociation cross section agrees 
with the result obtained independently by Brambilla et 
al. in their effective field theory approach [l9|, [2(| in the 
limit discussed in EDl . 



FIG. 3. (Color online) Top: Total decay width r to t (solid) and 
the partial widths for collisional damping Tdamp (dot-dashed) 
and gluodissociation Tdiss (dashed) for the T(1S) state. Bot- 
tom: Ratios of the partial widths rdi ss /r t ot (dashed) and 
Tdamp/rtot (dot-dashed) to the total width, and ratio of the 
partial widths Tdamp/rdiss (dotted). While collisional damp- 
ing is the dominant process in the QGP, gluodissociation can 
not be neglected for the T(IS') state. 



To obtain the mean gluodissociation cross section, 
we average our calculated gluodissociation cross sections 
over the Bose-Einstein distribution function of gluons at 
temperature T, thus assuming that the medium is ther- 
malized, although the heavy bb is not (see Fig. [2] for the 
gluon distribution): 
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where g e i — 16 is the number of gluonic degrees of free- 
dom. Taking Tdamp from the previous section together 
with the resulting width from gluodissociation yields the 
total decay width in the QGP, 
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Ttot as well as the partial widths Tdamp and Tdiss are 
displayed in Fig. [3] The ground-state width from col- 
lisional damping is seen to be about twice as large as 
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gluon-induced dissociation in the temperature range 200- 
400 MeV, such that both processes need to be con- 
sidered when calculating the total width in the quark- 
gluon plasma. Whereas damping increases monotonically 
with temperature, gluodissociation reaches a maximum, 
and decreases again at very high temperatures beyond 
600 MeV due to the diminishing overlap of the thermal 
gluon distribution and the gluodissociation cross section 
at large values of T. 



IV. TIME EVOLUTION OF THE FIREBALL 
AND DECAY CASCADE 

The density distribution of the lead ions is modeled 
by a Woods-Saxon potential with radius R = 6.62 fm 
and diffuseness a — 0.546 fm (26[. The number N bb - of 
produced 66-pairs at the point (x,y) in the transverse 
plain and impact parameter b is then proportional to the 
number of binary collisions JV C oi and nuclear overlap T AA , 
N bb (b,x,y) cx N con (b,x,y) cx T AA (b,x,y). The initial 
temperature is parametrized depending on the number 
of collisions, and Bjorken scaling is used for the time 
evolution [2 



T(b,t,x,y)=T c 



T AA (b,x,y) ( t, 
W0, 0,0) 



t 



1/3 



(7) 



where £qgp is the maximum lifetime of the quark-gluon 
plasma. We define a preliminary suppression factor R V AA , 
which accounts only for the bb suppression due to the 
three processes Debye screening, collisional damping and 
gluodissociation, 
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The numerator of eq. ([5]) is proportional to the number 
of bb bound states which have survived from their for- 
mation time tF until the fireball has cooled below the 
critical temperature T c , where the decay width Ttot is 
set to vanish. The integrand in the numerator of R AA , 



T AA {b,x,y)e 



J dtr tot (b,t,x,y) 
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is displayed in Fig. |4]for the T(IS') and T(2S) states for 
a central (6 = fm) and a peripheral collision (b = 10 
fm). The total density of the PbPb-system in the mo- 
ment of the collision where the nuclei pass through each 
other is also displayed. Clearly the T(25) is suppressed 
much more efficiently than the more stable T(15). Also 
one should note the action of Debye screening which for- 
bids the formation of bound bb states at sufficiently high 
temperatures and thus changes the shape of the surface 
from bell-shape (peripheral) into volcano-like (central). 
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FIG. 4. (Color online) Scaled populations (explanation in the 
text) of T(1S) (top) and T(2S) (middle) which remain after 
the fireball has cooled, projected on the transverse plane, for 
b — fm (left) and b = 10 fm (right) at £qgp = 6 fm/c 
and tF = 0.1 fm/c. The corresponding maximum density of 
both Pb-nuclei during the collision is also displayed (bottom) . 
Clearly the T(2S) is suppressed much more efficiently by the 
dissociation processes in the QGP. The suppression is stronger 
in the central regions of the collision where the temperature 
is higher. 



Results for the preliminary suppression factor of all 
five states for formation time tp = 0.1 fm/c and £qgp = 
6 fm/c are presented in Fig. [5] For the excited states 
Xt(2P) and T(3S) and higher excitations there exist no 
bound states for T > T c so that R A r A is equal for all these 
states. 

Now that we have calculated the suppression during 
the evolution of the fireball we have to consider the feed- 
down of the remaining bb population to calculate the frac- 
tion of decays into dimuon pairs, T(nS) —> . Fig. [5] 
displays the decays within the bb family and into dimuon 
pairs that are measured. Considering first the processes 
inside the fireball and then performing the decay cascade 
as a subsequent step is justified by the very different time 
scales involved. At the LHC the fireball has cooled within 
less than < 10 fm/c, while the subsequent decays take 
place on time scales ~ 10 3 fm/c. 

Let us denote bb states by / = (nl) and (C/j) (/ < 
J) the branching ratio of state J to decay into state / 
including all indirect decays with intermediate bb states. 
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FIG. 5. (Color online) Preliminary suppression factors 
H^ r ?(nl) from eq. © as functions of centrality for the dif- 
ferent bottomium states T(15*) (solid), T(25) (dash-dotted) 
Xb(lP) (dashed), and higher excited states (dotted) for the 
formation time tF = 0.1 fm/c and QGP lifetime £qgp = 6 
fm/c. 




FIG. 6. (Color online) Branching ratios for decays within 
the bottomium family T(nS) and \b{nP) and into /x ± -pairs 
according to [2S| . 



FIG. 7. (Color online) Suppression factor Raa for the T(1S) 
ground state calculated for 2.76 TeV PbPb-collisions from 
screening, collisional damping, gluodissociation and reduced 
feed-down using three QGP lifetimes tQGP = 4, 6, 8 fm/c 
(dotted, solid and dashed line respectively) for the centrality 
bins 50-100%, 40-50%, 30-40%, 20-30%, 10-20%, 5-10%, 
0-5%. The dash-dotted upper line is the preliminary sup- 
pression factor Ry^X (IS) (£qgp = 6 fm/c) without reduced 
feed-down. The corresponding CMS data [12] are in good 
agreement with the model results for the Y(IS') state. 



where B(nS — > ij^) is the corresponding branching ratio. 

We take N^ ± (nS) from the 2012 CMS data [3 and 
consider that 27.1% and 10.5% of the T(15) population 
comes from Xfe(l-P) and Xb(2P) decays, respectively [29j. 
Since these CDF results are obtained from pp collisions 
at 1.8 TeV with a transverse momentum cut p^ > 8.0 
GeV/c, it would be desirable to confirm the T(IS') pop- 
ulations from Xb decays in new pp measurements at 2.76 
TeV, which are not yet available. 

The initial populations are then obtained in units of 
B(nS -> ^NfpilS) as N^IS) = 16.3, N 1 (1P) = 
52.0, JV i (2S') = 20.7, N l {2P) = 39.4, ^(35) = 18.8. 
The final suppression factor is now simply calculated as 



R,AA{nS) 



NL(nS)/NUnS) 



The initial and final bb numbers of state /, N l (I) and 
N*(I) in pp and PbPb collisions are then connected by 

i<j 

< b PbW = £c„ivv)i^)- (io) 
i<j 

Further we define the number of T(nS) states that decay 
into dimuon pairs 

Nl ± (nS) = B(nS -> M ± )^ bPb (n5), (11) 



Raa{tiS) = B(nS — > (i = . 

Z2nS<J C U N ( J > 

(12) 



V. RESULTS 

We present the results for screening and collisional 
damping derived from the solutions of the Schrddinger 
equation with the potential eq. @, and the widths for 
gluodissociation as derived from eq. ([5]) . The total decay 
widths r to t are then inserted into a dynamic calculation 
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FIG. 8. (Color online) Suppression factors Raa for the T(25), Y(3S) states (top left and bottom left) and the double ratios 
(nS/lS)pbPb/(roS/lS) pp for n — 2,3 (top right and bottom right) calculated for 2.76 TeV PbPb-collisions from screening, 
collisional damping, gluodissociation and feed-down using three QGP lifetimes iQGP — 4, 6, 8 fm/c (dotted, solid and dashed 
line respectively) for the centrality bins 50-100%, 40-50%, 30-40%, 20-30%, 10-20%, 5-10%, 0-5% (left to right). The 
corresponding CMS results for the T(2S) state show significantly more suppression, in particular, in the peripheral region. 



for the fireball evolution to calculate preliminary suppres- 
sion factors, eq. (jHJ. 

Subsequently, the bottomium states pass through a de- 
cay cascade (see Fig. |6|) so that the higher excited states 
feed the lower lying states to yield the final suppression 
factor eq. (fT2l) . 

Our results for the suppression of the T(15) state in 
PbPb relative to pp are shown in Fig.[7]for three different 
QGP lifetimes iQGP = 4, 6, 8 fm/c as functions of cen- 
trality (number of participants). When comparing with 
our result from the preliminary suppression factor (up- 



TABLE I. Calculated minimum bias results for different tqap 
and tp = 0.1 fm/c compared to the CMS results [12J with sta- 
tistical and systematic error bars, respectively. The Raa(IS) 
is in good agreement with experiment, but the results for the 
excited states allow for additional suppression mechanisms. 



iQGP (fm/c) 


4 


6 


8 


CMS data [12] 


Raa{1S) 


0.52 


0.46 


0.42 


0.56 


±0.08 


± 


0.07 


Raa{2S) 


0.33 


0.28 


0.25 


0.12 


±0.04 


± 


0.02 


Raa(3S) 


0.28 


0.24 


0.22 


0.03 


±0.04 


± 


0.01 


(2S/lS) PbPb 
(2S/1S) PP 


0.63 


0.61 


0.61 


0.21 


±0.07 


± 


0.02 


(3S/lS) PbPb 
(3S/lS)p P 


0.53 


0.52 


0.52 


0.06 


±0.06 


± 


0.06 



per dotted step function), it is evident that the consider- 
ation of the feed-down cascade is essential for modeling 
the suppression. 

The CMS data point [12] at 40-50% centrality violates 
the monotonic increase of the suppression with centrality, 
but is consistent with the other points within statistical 
and systematic error bars. Hence, the calculated sup- 
pression is in very good agreement with the CMS data 
for the T(1S) ground state. This is also true for mini- 
mum bias (centrality integrated) results, which are shown 
in Table |U 

Our results for the suppression of the T(25) and T(3S) 
states in PbPb relative to pp are shown in Fig. [5] (left 
column) for three different QGP lifetimes tQGP = 4, 6, 8 
fm/c as functions of centrality. The double ratios with re- 
spect to the T(15) state in PbPb and pp are displayed in 
the right column of Fig. [SJ with CMS data [3 included 
for the T(25) state. The suppression found experimen- 
tally for the T(2S) state is much more pronounced than 
in the calculation, in particular, for the three more pe- 
ripheral data points. 

It appears to be very difficult for theoretical models to 
obtain such a huge suppression of the T(2S) state in pe- 
ripheral collisions, and indeed, other approaches such as 
[30l - [32| also find that the T(2S) suppression factor rises 
towards 1 for peripheral collisions. As a consequence of 
the disagreement with the centrality-dependent data, our 
minimum-bias results of Table U also disagree substan- 
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tially for the T(2S) and T(3S) states. 

The reason for the disagreement will probably be 
cleared up once more precise pp reference data at 2.76 
TeV become available in the future. It is, however, also 
conceivable that additional suppression mechanisms not 
considered in this work play a role for the T(2S) and 
T(35) states. 



VI. CONCLUSION 

We have formulated a three-step model for the suppres- 
sion of the bottomium states T(nS) in the quark- gluon 
plasma that is formed in PbPb collisions at LHC energies. 
Due to its stability against screening up to very high tem- 
peratures, the T(IS') state is a particularly suitable probe 
for the relevance of gluodissociation, collisional damping, 
and reduced feed-down. 

We find that gluodissociation of the Y(LS') state is size- 
able [13j due to the strong overlap of the T(15) gluodisso- 
ciation cross section with the thermal gluon distribution. 
In the temperature region 200-400 MeV, both gluodis- 
sociation and collisional damping are found to be impor- 
tant. 

The observed suppression factor Raa(1S) — 0.56 in 
minimum-bias PbPb collisions [l2[ is essentially due to 
gluodissociation and damping of the T(15) state, and to 
the melting and dissociation of the excited states: The 
excited states in particular, the Xb{nP) states par- 
tially feed the T(1S) state in pp, pp and e + e~ collisions, 
and their melting and dissociation in the quark-gluon 
plasma substantially reduces the feed-down in PbPb col- 
lisions at LHC energies. 

The calculated T(15) suppression factor as function of 
the collision centrality is indeed in very good agreement 



with the CMS data if the modification of the feed-down 
cascade in PbPb as compared to pp is taken into account. 

Different from the T(15) ground state, the excited 
states - and in particular, the T(25) and T(3S) states 
that are observed in the CMS experiment - are already 
suppressed through screening to a much larger extent 
than the ground state, so that the contributions from 
damping and gluodissociation are less important here. 
The dissolution of the excited states in the quark-gluon 
plasma causes the substantial feed-down reduction that 
is one of the three main reasons for the ground-state sup- 
pression. 

From our calculations it appears that there may be ad- 
ditional causes for the suppression of the excited states, 
such as cold nuclear matter (CNM) effects - although 
these should essentially cancel out in the double ratios 
that are shown in Fig. [8] It is conceivable that CNM- 
effects will be constrained in forthcoming pPb measure- 
ments at the LHC. Compared to the present CMS ex- 
perimental results for the suppression of the T(25) and 
T(35) states in PbPb [HI, our calculated Raa values 
are substantially too large, in particular, in peripheral 
collisions. 

Apart from these detailed deficiencies, it appears that 
the overall good description of Upsilon suppression in 
PbPb collisions at LHC energies in the present formu- 
lation contributes to the available indirect evidences for 
quark-gluon plasma formation in heavy-ion collisions at 
high relativistic energies. 
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